Scaling of Crack Surfaces and Implications on Fracture Mechanics 
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The scaling laws describing the roughness development of crack surfaces are incorporated into 
the Griffith criterion. We show that, in the case of a Family- Vicsek scaling, the energy balance leads 
to a purely elastic brittle behavior. On the contrary, it appears that an anomalous scaling reflects 
a TJ-curve behavior associated to a size effect of the critical resistance to crack growth in agreement 
with the fracture process of heterogeneous brittle materials exhibiting a microcracking damage. 
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Fractography has always been a useful tool for the un- 
derstanding of fracture mechanisms in heterogeneous ma- 
terials. It is well established that the fracture properties 
of brittle heterogeneous materials such as rock, concrete, 
tough ceramics, wood and various composites strongly 
depend on the microstructure and on the damage pro- 
cess of the material. Both these parameters are also at 
the source of the roughness of crack surfaces. Indeed, 
microstructural heterogeneities acting on the crack front 
as obstacles, and the existence of a large fracture process 
zone with cracking damage where elastic interactions oc- 
cur between cracks, may have a strong influence on the 
local deviations of the main crack |jj-p|]. Thus, the rough- 
ness of fracture surfaces can be considered as an inheri- 
tance of the heterogeneous character and of the damage 
process of the material. Hence, it is naturally tempting 
to correlate the fracture surface morphology to macro- 
scopic mechanical properties such as fracture energy or 
toughness. 

The studies of quantitative fractography are nowadays 
a very active field of research. Many experiments |7|-|lC|l 
on materials as different as ductile aluminium alloys to 
brittle materials like rock or wood have shown that the 
topography of fracture surfaces is self-affine fTl}] . How- 
ever, in most cases and in spite of strong differences in the 
materials, a roughness index £/ oc = 0.8 — 0.9 (called local 
roughness exponent in what follows) has been reported in 
all cases. This robustness of the results seems to support 
the idea suggested by Bouchaud et al. || that the local 
roughness exponent might have a universal value, i.e., 
independent of the fracture mode and of the material. 

Nevertheless, if the local roughness exponent is univer- 
sal, the range of lengthscales within which this scaling 
domain is observed strongly depends on the material mi- 
crostructure. Thus, recent studies focussed on the com- 
plete description (3D) of the crack morphology of granite 
H and wood H have shown that the scaling laws govern- 
ing the crack developments in longitudinal and transverse 
directions are different and material dependent. Let us 
consider the development of a fracture surface from a 
border where a straight notch of length L imposes a zero 
roughness. The mean plane of the crack surface is de- 
fined as (x, y) where the x axis is perpendicular to the 
direction of crack propagation and the y axis is parallel 



to this propagation direction. It has been found that the 
height fluctuations Ah of fracture surfaces of two het- 
erogeneous brittle materials (granite j§| and wood 
estimated over a window of size I along the x axis and at 
a distance y from the initial notch exhibited intrinsically 
anomalous scaling properties obtained in some models of 
nonequilibrium kinetic roughening mM : 



Ah(l,y)~A 
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where £(y) — By 1 /* depends on the distance to the initial 
notch y and corresponds to the crossover length along the 
x axis below which the surface is self-affine with the local 
roughness exponent Q oc . Along the y axis, the rough- 
ness develops according to two different regimes : for 
large lengthscales, i.e. I 3> £(y), the roughness grows 



as Ah(l, y) 



If 



where £ is called the global rough- 



ness exponent while for small lengthscales, i.e. I <C £(y), 
the roughness growth is characterized by the exponent 
(C ~ Cioc)/z. The global roughness exponent £, the expo- 
nent z (called dynamic exponent) and the prefactors A 
and B are material dependent. Thus, despite the univer- 
sal distribution of the height fluctuations in the trans- 
verse direction of crack propagation, the roughening in 
longitudinal direction is material dependent and has an 
influence on the roughness magnitude in the transverse 
direction. All lengths are dimensionless. Real lengths are 
divided by a unity length I*. 

One main consequence is that, when the global satura- 
tion occurs, i.e. far from the notch (for y 3> y sa ti where 
y sa t = (L/B) z ), the magnitude of the roughness is not 
just a function of the window size but also of the system 
size L : Ah(l,y 3> y sa t) — Al^ loc L^~^ loe . This unconven- 
tional dependence of the stationnary local fluctuations on 
the system size is distinctly different from what happens 
in the Family- Vicsek scaling |l2j which is defined as : 
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The Family- Vicsek scaling can be seen as a particular 
case of anomalous scaling where: C — C/oc- 

In this study, we propose to discuss the link between 
the roughening of crack surfaces and the fracture process. 
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On the basis of the Griffith criterion, we show that a frac- 
ture surface exhibiting an anomalous scaling reflects fun- 
damental mechanical behavior of quasi-brittle materials: 
i2-curve behavior (i.e. an evolution of the resistance to 
crack growth as a function of the crack length increment 
||) and size effect of the elastic energy release rate. Fig. 1 
shows an illustration of a i?-curve behavior observed for 
a wood crack propagation (mode I) using a TDCB tech- 
nique ||[l3| . The energy release rate G (proportional to 
the square of the toughness) evolves very significantly at 
the onset of crack propagation and becomes independent 
of the crack length after a characteristic propagation dis- 
tance. In Fig. 2, an example of size effect on this critical 
energy release rate at saturation is presented [p^|] . 

Several connections between the fractal dimension of 
the crack surfaces and the fracture energy or toughness 
have been proposed [pp|Jl9|-p^| . Nevertheless, most of 
the proposed analytical models |l9|-|2~4|] show weak con- 
nections with the fracture behavior of brittle heteroge- 
neous materials and especially with the phenomenolog- 
ical i?-curve behavior and size effect. However, these 
models have been built on the basis of a 2D self-afhnc 
crack surface disregarding the anisotropy of the scaling 
laws governing the crack developments in longitudinal 
and transverse directions. In this case, either a ii-curve 
is obtained without any size effect or a size effect with- 
out i?-curve. Let us consider a semi-infinite linear elas- 
tic specimen of thickness L containing an initial crack 
at position Aa submitted to an uniaxial stable and slow 
tension (mode I). In the classical Griffith approach, the 
fracture criterion is estimated by balancing the elastic 
energy released at the macroscale during an infinitesimal 
crack propagation and the energy required to create the 
corresponding free surfaces at the microscale. According 
to linear elastic fracture mechanics (LEFM), the rate of 
energy dissipation in the structure as a whole must be 
defined with respect to the projected crack surface A p 
(at macroscale, crack roughness can be ignored). On the 
contrary, the estimate of the energy required for crack 
propagation at the microscale needs to take into account 
the real crack surface A r . Introducing the elastic energy 
release rate G at the macroscale, the fracture criterion 
becomes for an infinitesimal crack advance : 



G SA P = 2 7 SA r 



(3) 



where SA r is the real area increment and 5A p its pro- 
jection on the fracture mean plane. The term 7 is 
the so-called specific surface energy that characterizes 
the resistance of material to cracking. The real area 
can be obtained from the surface description during the 
crack growth da from position Aa of the crack front : 

6A r ~ f&a +5a i'iv) dy, where tp(y) is the length of a pro- 
file at position y of the crack roughness. For a specimen of 
width L, a first order approximation leads to the follow- 
ing expansion of the fracture criterion: G = 2jip(y)/ L. 



According to (|]), the roughness Ah(l,y) can be seen as 
a description of the fracture profile at position y. The 
length i\) can be estimated by covering the profile path 
with segments of length 5 the horizontal projection on the 
x axis is l and the vertical one on the z axis is on average 
[seeEq.(0)] : A^y)^^ l o 0oc . The length l a is chosen as 
the lower cutoff of the fractal range ( i.e. the character- 
istic size of the smaller micro-structural element relevant 
for the fracture process). Moreover, the crack profile is 
a univalued function (observed fracture surfaces do not 
show overhangs) and hence the number of segments 8 
needed to cover the curve is given by the ratio (L/l ). As 
a consequence the real length ip = (L/l )5 of the crack 
profile appears dependent on the distance y to the initial 
notch and can be expressed as : 



(1-Cloc 
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In Eq.(^), the terms between brackets correspond to the 
tangent of the typical angle between the segment 5 and 
the x axis. Thus, as expected intuitively, the length of 
a self-affine curve is larger than its projected length al- 
though it is proportionnal to it. On the other hand, ac- 
cording to (||) when the global saturation state of the 
roughness is reached (i.e. for y 3> y sa t), the length of 
the crack profile saturates and becomes a nonlinear func- 
tion of the specimen width L. Note that the length of 
the self-affine crack profile [Eq.(S)] is a finite quantity 
in contrast with the expressions proposed by Carpinteri 
ifjff , Borodich (§3] and Bazant [Q. In these models, the 
actual length of a fractal curve is considered as infinite 
and this leads the use of an unconventional definition of 
fracture energy (fractal fracture energy) . As we will show 
in the following, the introduction of a lower cutoff for the 
surface description allows to obtain a fracture energy in 
agreement with the classical dimensions of LEFM. 

In the zone where the roughness grows (i.e. for Aa <C 
Aa SQ t with Aa sat = y sa t), the fracture equilibrium (||) 
leads to an energy release rate function of Aa : 



G R (Aa < Aa sat ) 




Aa 2(C-£ l0 c)/* 



(5) 



In Eq.(|^), the subscript R emphasizes that the behav- 
ior of the resistance to fracture growth is similar to a 
i?-curve Q. When the crack increment is large, i.e. for 
Aa 3> Aa sat which corresponds to the saturation state of 
the roughness, the resistance to fracture growth becomes 
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G RC (Aa » Aa sat ) ~ 2 7 ^ 1 + (^TT^rJ ^^'^ 

(6) 

In the zone where the roughness saturates, the length of 
the crack profile remains constant which induces a resis- 
tance to crack growth independent of the crack length 
increment Aa (^J). The subscript C in Grc empha- 
sizes that the resistance to crack growth has reached an 
asymptotic or critical value. This post i?-curve regime 
at constant resistance to crack growth obtained here in 
an analytical form [Eq.(^)] is in agreement with the one 
generally assumed in the phenomenological approaches 
(for instance p5[). In Fig. 1, Eq.(||) and Eq.(^|) are used 
to fit the data. Eq. (||) provides a good description of the 
growth of the energy release rate. This fit needs to use 
an overestimated ratio A/l l ~^ lo ° . Note that, a previous 
approach obtained from the energy released by a frac- 
tal pattern of microcracks ahead of the main crack of a 
quasi-brittle material [B3[ have shown equally a i?-curve 
behavior. On the other hand, it has been shown that the 
roughness of a crack induces a decrease of the stress in- 
tensity factor and this shielding process is also known 
to be a source of the phenomenological i?-curves. These 
latter points emphasize that the connection between the 
anomalous roughening of fracture surfaces and fracture 
mechanics seems to reflect the particular fracture behav- 
ior of brittle materials showing toughening mechanisms. 

The main consequence of the anomalous scaling is the 
dependence of the critical resistance Grc to crack growth 
on the specimen size L [see Eq.(||)]. Size effect phe- 
nomena is one of the main characteristic of the fracture 
behavior of heterogeneous brittle materials such as con- 
crete, rocks or wood |2j| . Figure || shows the theoreti- 
cal evolution of critical resistance to crack growth as a 
function of the system size (or characteristic size) L ob- 
tained for arbitrary values. It appears that the size effect 
shows two asymptotic behaviors respectively Grc = 2 7 
and Gb.c ~ L^~^ oc . From Eq.(||), one can define the 
crossover length L c = (^ O 1_ ^ oc /A) 1 /('*~'* , ° c \ p or sma n 
system sizes L <C L c , which correspond theoretically to 
"shallow" surfaces, i.e. surfaces that have a mean local 
angle of the crack profile smaller than 45 degrees, there 
is no size effect : Grc — 2 7 . The roughness of crack 
surfaces being negligible, real crack surfaces are weakly 
different from the projected one and so, classical results 
of purely elastic fracture mechanics are recovered. For 
large system sizes L 3> L c , the critical resistance evolves 
as a power law : Grc ~ L^~^ ,oc . These fracture surfaces 
are theoretically "spiky" (i.e. the average local angle of 
the crack profile is greater than 45 degrees) and in this 
case, the critical energy release rates are greater than the 
specific surface energy 2-f. Note that the latter power law 
behavior is in good agreement with the observed size ef- 
fect in wood as shown in Fig. 2. In fact, the typical 



local angle of wood crack surfaces is smaller than 45 de- 
grees but the crossover length L c , calculated with the 
ratio A/l 1 ~^ l ° c obtained from the i?-curve fit (Fig. 1), is 
smaller than the sizes L of tested specimens. 

If the development of the roughness crack is driven by 
a Family- Vicsek scaling instead of an anomalous scaling, 
the crack profile length ip is independent of the crack in- 
crement Aa: -0(y) ~ 1 + (A / / O 1-Cioc ) 2 L. Hence, the 
resistance to fracture growth reduces to : 

G c (Aa)^2 1] j^(^f (7) 

Thus, for a Family- Vicsek scaling, the resistance to frac- 
ture is independent of the crack position Aa and of the 
specimen size L. This is consistent with the behavior of 
purely elastic brittle materials. Likewise, if the global 
roughness exponent tends to the local one (£ — * (i oc ), 
the i?-curve behavior (|B|) and the size effect vanish 
progressively and the fracture behavior is close to purely 
elastic brittle materials [Eq.(Q)]. 

In conclusion, on the basis of complete descriptions 
(3D) of fracture surface morphologies given by Family- 
Vicsek and anomalous scaling laws, we have shown that 
the connections between the roughening of fracture sur- 
faces and fracture mechanics are important. A Family- 
Vicsek scaling of the fracture roughness reflects a purely 
clastic brittle fracture behavior, while the fracture be- 
havior obtained on the basis of an anomalous scaling ac- 
counts for a i?-curve behavior and size effect of the criti- 
cal resistance to crack growth. Anomalous scaling reflects 
the experimental behavior of brittle heterogeneous mate- 
rials showing toughening mechanisms. On going experi- 
mental studies attempt to link systematically mechanical 
behaviors and fracture roughening for different materials. 

We are particularly greatful to J.R. Rice for fruitful dis- 
cussions and enlighting remarks. This work was partly 
supported by the CNRS-ECODEV project. 
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FIG. 1. i?-curve behavior observed for a TDCB fracture experiment in wood. For this Norway Spruce specimen of charac- 
teristic size L — 30mm, data are fitted from Eq.(|^) and Eq.Q with the scaling exponent (£ — C,i oc )/z = 0.32 which is consistent 
with the expected result for this specimen : (£ — C,i oc )/z — 0.36 ± 0.09. 

FIG. 2. Size effect on the critical energy release rate Grc (Norway Spruce specimens, 11.3 < L < 60mm). The behavior 
derived from the anomalous scaling (solid line): Grc ~ L^ 1 ' 100 [Eq.Q] is characterized by the exponent 0.69 which is in good 
agreement with the expected result for Spruce : ( — £( oc = 0.73 ± 0.17. 



FIG. 3. Theoretical size effect on the critical resistance to crack growth [Eq|6J obtained from an anomalous scaling and for 
the arbitrary scaling exponents : £ = 1.3, 0oc = 0.8, A — 0.1 and l = 1. 
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